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Abstract 



Recently, it has been shown that properties of excitable media stirred by two- 

<^ • dimensional chaotic flows can be properly studied in a one-dimensional framework 

CD . [8,9], describing the transverse profile of the filament-like structures observed in the 

\ system. Here, we perform a bifurcation analysis of this one-dimensional approxima- 

■ tion as a function of the Damkdhler number, the ratio between the chemical and 

. the strain rates. Different branches of stable solutions are calculated, and a Hopf 

' bifurcation, leading to an oscillating filament, identified. 

^. 
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1 Introduction 



Chemically reacting substances in fluid flows are complex systems important 
both from a fundamental point of view and for its relevance in industrial and 
environmental contexts [1,2]. Typical incompressible chaotic flows stretch fluid 
parcels along particular directions, with the consequent contraction along the 
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others. Thus, passive markers draw filamental or lamehar structures that have 
been investigated both from theory as from experiment [3,4]. Chemical reac- 
tants are also stretched in this way, producing a great increase in the surface 
of contact between different species with deep effects on the global chemical 
kinetics [5]. In addition to the strictly chemical interactions, biological interac- 
tions in ecosystems (predation, grazing, consumption, competition, ...) can also 
be described in the same framework, so that stability and dynamics of aquatic 
ecosystems are also strongly influenced by this filamentation process[6]. 

In this Paper we analyze a simplified model describing the transverse chemical 
structure of these filaments, for the case in which the chemical or biological ac- 
tivity is of the excitable type. This includes in particular well known laboratory 
reactions, such as the Belousov-Zhabotinsky reaction, but also reactions of en- 
vironmental importance, such as phytoplankton-zooplankton competition [7]. 
Such kind of dynamics in two-dimensional chaotic fiows was analyzed in [8,9]. 
There it was found that some of the dynamic regimes and transitions be- 
tween them can be understood in terms of a simple one-dimensional model, 
of the kind already considered in [10,11,12,13] that focusses in the transverse 
structure of the filaments. It was mentioned in [9] that, in addition to the 
bifurcations considered there in detail, there was a range of parameters for 
which a complex coexistence of solutions and bifurcations occurred, that may 
be of relevance to understand qualitative behaviors of the two-dimensional 
chaotic flows. The present paper focusses in that regime. We use the FitzHug- 
Nagumo model as a prototypical excitable dynamics, but we expect that the 
same qualitative results will be also obtained for other chemical or biological 
excitable schemes. Similar behavior has been observed in an excitable model 
of oceanic plankton population in [14] . 



2 The one-dimensional filament model 



Given a set of kinetic laws describing the time evolution of a number of chem- 
ical or biological concentrations {Cj(x, jv in a homogeneous system, 

^ = A;F,(Ci,...,C^) (1) 



{k is a global reaction rate), a general partial differential equation model de- 
scribing their spatiotemporal evolution under the simultaneous additional ef- 
fects of advection in an incompressible velocity field v(x, t), and diffusion with 
diffusion coefficients {Di}, is the one given by the following advection-reaction- 
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diffusion equations: 



dCi{^,t) 
dt 



+ v(x,t)-Va(x,i) = A;F,(Ci,...,Cjv) + AV2C,(x,i) . (2) 



In tfie following we restrict for simplicity to the case of equal diffusion coef- 
ficients Di = D, In the immediate vicinity of filaments or lamelae, one 
can approximate the fiow by its linearization around a point co-moving with 
a fluid element. In areas where the advection dynamics is hyperbolic this lin- 
earization identifies principal directions {xj} along which the corresponding 
velocity components read Vj — ^jXj. Positive and negative values of the strain 
rate identify expanding and contracting directions, respectively. Along the 
expanding directions, diffusion and advection cooperate and the chemical con- 
centrations are fastly homogenized, whereas strong gradients build-up in the 
contracting directions. This is the origin of filamental (one expanding direc- 
tion) or lamellar (two expanding directions) structures. Gradients in Eq. (2) 
can be safely neglected except along the contracting directions. We consider 
here the case of filaments in two-dimensional fiows, or lamellae in three dimen- 
sions, so that there is only one contracting direction, of strain rate Aj = —A, 
with A > 0. In this case, (2) reduces to an effective one-dimensional model 
for the transverse profile of the chemical distributions. By measuring time in 
units of A~^, and space in units of JD/\ we arrive at: 



We have introduced a Damkohler number Da = A;/ A as the ratio between the 
chemical and the strain rate. 

Several limitations are inherent to (3). It is a strictly local model valid close 
to (moving) hyperbolic trajectories. The linearization assumption will only be 
justified if the full filament width is smaller than any structure in the velocity 
field. This implies in particular a very small diffusion coefficient so that we are 
deeply in the Batchelor regime. Additional qualitative changes of behavior, not 
contained in (3), are expected when increasing the diffusion coefficient in (2). 
Curvature effects or, more importantly, any interaction between filaments, are 
completely neglected. Finally, the use of a constant A neglects any fiuctuations 
of the stretching rate. 

The FitzHugh-Nagumo (FN) model consists in a dynamics of the type (1) for 
two interacting species, of concentrations Ci and C2, and reaction terms: 



dt 




(3) 




(4) 
(5) 



F2 = e(Ci-7C2). 
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Fig. 1. One- hump solutions for the two chemical components: Ci and C2 (multiplied 
by a factor 5); a) Da = 25; b) Da = 55. 

The FN model behaves excitably when e <S 1 so that there is a separation 
between the fast evolution of the active component or activator, Ci, and the 
slow evolution of C2, the passive one or inhibitor. We concentrate in the pa- 
rameter values e = 10~^, 7 = 3.0, and a — 0.25, for which robust excitable 
behavior is obtained. 



3 Filament solutions 

The unexcitcd solution, Ci = C2 — 0, is a linearly stable exact solution of 
system (3)- (5). In addition, the most notable solution is a pulse like steady 
solution, centered on a; = 0, in which the activator is fully excited in a central 
region, whereas the inhibitor remains at reduced concentrations. This chemical 
distribution can be thought as the transverse cut of one of the excited filaments 
seen in two-dimensional excitable fluid flows. Pulse solutions of this type are 
shown in Fig. 1. 



The width, Wg, of this solution is determined by the competition between the 
tendency to expand, due to the joint effects of the excitable chemistry and 
diffusion, and the compression by the flow, and can be estimated [8,9] to be, 
for small e, ~ (1 — 2a)y^Da/2. It is very similar to what is obtained in one- 
component bistable chemical models [12], so that the presence of the inhibitor 
plays here only a secondary role. 

As in bistable chemical models, this stable pulse solution coexists in a range 
of parameters with a unstable pulse, of width given at small e approximately 
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Fig. 2. Two types of pulses for Da = 100: a) The two-hump symmetric solution; b) 

the asymmetric one. 

by [9] Wu ~ 2/vaDa. When the value of Da decreases (meaning that the 
chemistry becomes slower, or the stretching faster) the two solutions approach, 
collide and disappear from phase space in a saddle- node bifurcation [9]. This 
happens at a value of Dac ~ 2y^2/a(l — 2a)~^. Thus, the excited pulse solution 
does not exists at small Da, and the existence of a critical Dac explains [9] 
the dynamical transition occurring in closed two-dimensional chaotic flows 
between a situation in which local perturbations of the unexcited state have 
limited impact on the system, and a state in which they give rise to a global 
excitation of the whole fluid. A related transition occurs in open flows [9]. 

At larger values of Da, the influence of the inhibitor becomes more notice- 
able and, as a result, steady filament solutions have a 'two-humped' shape 
(Fig. 2a). It was noted in [9] that this kind of 'double-filament', also seen in 
two-dimensional flow simulations, can be regarded as a bound state of two 
asymmetric counter-propagating excitable pulses, which are also steady solu- 
tions of model (3)-(5) coexisting with the symmetric ones at large Da (Fig. 
2b). Of course, for each value of Da there are two asymmetric solutions of this 
kind, each one being the specular image of the other. 



4 Bifurcation behavior at intermediate Da 



It was noted in [9] that a complex bifurcation scenario occurred between the 
small and the large values of Da for which the solutions in Figs. (1) and (2) 
were respectively obtained. It happens that these three solution branches are 
not directly connected, at least for the values of e and 7 considered here. We 
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Fig. 3. Bifurcation diagram for the system given by Eq. (3). We plot the height of 
stable solutions in terms of Da. Solid line corresponds to decreasing the value of 
Da from an initial state of the asymmetric type, like that in Fig. 2b), at large Da. 
Dashed line is obtained in the same way but beginning from a symmetric two-hump 
pulse (Fig. 2a). Dashcd-dotted line is obtained increasing the value of Da from a 
symmetric one-hump pulse (Fig. la)). 



have followed the three branches of steady stable pulse solutions by starting 
with well developed solutions at large and small Da, and then performing slow 
changes in Da. In this way we can only follow stable solutions. A more complete 
characterization would require also the continuation of unstable branches, that 
will be performed elsewhere. We monitor the height of the activator at x — 0, 
and plot it in Fig. 3. 



When increasing the value of Da on the symmetric one-humped solution 
branch, the inhibitor value increases in the center, so that for Da > 45, 
the central concentration of the activator becomes a shallow minimum (Fig. 
lb). So, above this parameter value, this branch can be called 'slightly two- 
humped', and it remains so until Da ~ 77.1. It is however clearly different 
from the 'strongly two-humped' solution branch to which Fig. 2a pertains. 
At Da ~ 77.1, a Hopf bifurcation occurs: the filament (both the Ci and the 
C2 concentrations) pulsates in height, width and shape. We show in Fig. (4) 
this behavior. For Da above this Hopf bifurcation, we plot in Fig. 3 the max- 
imum and minimum central height values attained during the oscillation. At 
Da 81.0 the width becomes too narrow at some moment of the oscillation 
and the filament collapses to the unexcited (Ci = C2 = 0) solution, so that 
the limit cycle solution disappears. This collapse probably reveals the colli- 
sion of the filament limit cycle with an unstable pulse solution that we have 
not characterized. The basin of attraction of the oscillating solution is rather 
narrow: at these values of Da it is easier to be attracted by the two-hump 
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Fig. 4. Spatio-temporal plot of Ci in an oscillating filament at Da = 79. 
symmetric steady solution branch or the unexcited state. 



On the other hand, if we start with a well developed two-hump filament 
such as the one in Fig. 2a at high Da, the two humps approach each other 
when decreasing Da. Before they fuse together, there is a discontinuous jump 
(at Da ^ 75.4) to the steady 'slightly two-humped filament' followed before. 
This probably reveals a saddle-node bifurcation in which the stable two-hump 
branch collides with an unstable one, that we have not followed. 

Finally, starting from the asymmetric filament of Fig. 2b, it approaches the 
X = axis when decreasing Da, becoming more and more symmetric. For 
the values of e and 7 used here, however, it does not join smoothly with the 
symmetric branch in a forward pitchfork bifurcation in which it would also 
join the other asymmetric filament of opposite symmetry. Rather, it performs 
a (small) discontinuous jump to the symmetric filament branch at Da ~ 70.48. 
This is probably the signature of a backward pitchfork bifurcation, or of some 
other complex coexistence with unstable solutions, that would also explain the 
observed range of bistability between the symmetric and asymmetric filament 
solutions. 
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5 Summary 



We have investigated several pulse-like solutions of a one-dimensional model 
intended to characterize the chemical structure transverse to filaments or lame- 
lae generated by chaotic flows in reacting systems. The most striking feature 
has been the finding of an oscillating solution, in which an activator-inhibitor 
competition is sustained forever by the compressing strain. Moreover, param- 
eter ranges have been found in which there is coexistence of up to three stable 
pulse solutions, in addition to the unexcited one. The transition behavior ob- 
served in two-dimensional open and closed chaotic flows at intermediate Da 
numbers [9] is probably related to this complex phase space. Further work is 
needed to fully characterize this relationship. 
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